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ABSTRACT 

Observations from Spitzer and ground-based infrared spectroscopy reveal significant diversity in the 
molecular emission from the inner few AU of T Tauri disks. We explore theoretically the possible 
origin of this diversity by expanding on our earlier thermal-chemical model of disk atmospheres. We 
consider how variations in grain settling, X-ray irradiation, accretion-related mechanical heating, and 
the oxygen-to-carbon ratio can affect the thermal and chemical properties of the atmosphere at 0.25- 
40 AU. We find that these model parameters can account for many properties of the detected molecular 
emission. The column density of the warm (200-2000K) molecular atmosphere is sensitive to grain 
settling and the efficiency of accretion-related heating, which may account, at least in part, for the 
large range in molecular emission fluxes that have been observed. The dependence of the atmospheric 
properties on the model parameters may also help to explain trends that have been reported in the 
literature between molecular emission strength and mid-infrared color, stellar accretion rate, and disk 
mass. We discuss whether some of the differences between our model results and the observations 
(e.g., for water) indicate a role for vertical transport and freeze-out in the disk midplane. We also 
discuss how planetesimal formation in the outer disk (beyond the snowline) may imprint a chemical 
signature on the inner few AU of the disk and speculate on possible observational tracers of this 
process. 

Subject headings: (stars:) circumstellar matter — (stars:) planetary systems: protoplanetary disks — 
(stars:) planetary systems: formation 
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1. INTRODUCTION 

Recent work with the Spitzer Space Telescope has 
found that emission from water and organic molecules 
is common in spectra of T Tauri stars (Carr & Najita 
2011, 2008; Salyk et al. 2011, 2008; Pontoppidan et al. 
2010; Pascucci et al. 2009), with the emission likely aris- 
ing from the inner few AU region of the circumstellar 
disk (Carr & Najita 2008; Salyk et al. 2008). These 
new molecular diagnostics complement the more famil- 
iar near-infrared molecular diagnostics (CO, water, OH) 
that probe disk radii < 1 AU (see Najita et al. 2007; Car- 
mona 2010 for recent reviews). Since T Tauri disks are 
optically thick in the continuum, these molecular emis- 
sion features presumably arise from a temperature inver- 
sion region in a warm disk atmosphere. The possibility of 
warm molecular emission from the atmospheres of inner 
disks is supported by thermal-chemical models of disk 
atmospheres (Glassgold & Najita 2001; Glassgold et al. 
2004; Kamp & DuUemond 2004; WiUacy & Woods 2009; 
Nomura & MiUar 2005; Gorti & Hollenbach 2008). 

One of the interesting properties of the mid-infrared 
molecular emission observed with Spitzer is that the 
fluxes and relative strengths of molecular features vary 
among the T Tauri stars (Carr & Najita 2011; Salyk 
et al. 20011; Pontoppidan et al. 2010; Pascucci et al. 
2009). The strength of the water emission varies by a 
factor of > 30 (e.g., Pontoppidan et al. 2010), and an 
equally large range is found for the strength of the HCN 
emission (e.g., Salyk et al. 2011). The water emission 
can be brighter or fainter than the HCN emission, and 



HCN is usually stronger than C2H2 among T Tauri stars 
but is sometimes weaker (Carr & Najita 2011; Salyk et 
al. 2011; Pontoppidan et al. 2010; Pascucci et al. 2009). 
What is the origin of these variations in flux and relative 
strength? 

Several interesting possibilities are suggested by trends 
in the demographics of molecular emission sources. 
While HCN is typically stronger than C2H2 in T Tauri 
stars, Pascucci et al. (2009) found C2H2 to be stronger 
than HCN in median spectra of brown dwarfs. Pascucci 
et al. suggested that this dichotomy arises because of the 
stronger UV fleld of T Tauri stars, which enhances the 
dissociation of N2 and the formation of HCN in those 
systems. 

Further possible trends in HCN emission strength were 
noted by Teske et al. (2010). They found that the HCN 
emission strength of T Tauri stars appears to increase 
with stellar accretion rate and with X-ray luminosity, 
suggesting that these parameters may play a role in 
the molecular emission fluxes. They suggested that the 
trend of HCN flux with stellar accretion rate might re- 
sult from enhanced accretion-related mechanical heating 
in the disk atmosphere, along the lines explored theo- 
retically by Glassgold et al. (2009). These authors had 
previously suggested that increased rates of accretion- 
related mechanical heating can increase the column den- 
sity of warm water in the T Tauri disk atmosphere, pro- 
ducing stronger water emission. A similar effect might 
also be expected for other molecular species. Consistent 
with this perspective, Salyk et al. (2011) have reported 
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a higher detection rate of molecular emission in systems 
with larger Ha equivalent widths (and tentatively also 
with stellar accretion rate), which they interpret as sug- 
gesting that accretional heating plays a role in enhancing 
the line emission. Salyk et al. also find that molecular 
emission is preferentially detected from disks that have 
experienced more grain settling, as indicated by their 
mid-infrared colors. 

As an alternate explanation for their results, Teske et 
al. also discussed the possibility that increased UV irra- 
diation produced at higher stellar accretion rates might 
enhance the HCN abundance, as discussed by Pascucci 
et al. (2009). Little interpretation was offered for the 
possible trend they found with Lx. other than that X- 
ray irradiation may help to synthesize HCN by raising 
the abundance of molecular ions and radicals. 

Carr & Najita (2011) also reported a possible relation 
between the HCN/H2O emission flux ratio and disk mass 
in their small sample of T Tauri star spectra. The sources 
with the largest HCN/H2O flux ratios are those with the 
largest disk masses. They note that such a trend might 
be expected as a consequence of planetesimal and pro- 
toplanet formation in the outer disk. As described by 
Ciesla & Cuzzi (2006 and references therein), because icy 
bodies ranging from centimeter to kilometer and larger 
size migrate at different rates, their migration can lead 
to enhancements or depletions of water vapor in the disk 
that may vary as a function of both disk radius and time. 
Since water is a signiflcant reservoir of oxygen, these ef- 
fects can lower or raise the gas-phase oxygen-to-carbon 
(0/C) ratio in the inner disk and possibly affect the HCN 
and H2O abundances in the manner observed. More 
generally, because it is difflcult to diagnose the forma- 
tion of kilometer-sized planetesimals, the building blocks 
of planets, by other means, it is interesting to explore 
whether their formation may induce this (or any other) 
chemical signature in the inner disk that may signal their 
formation. 

Here we explore the above ideas by expanding on 
our earlier thermal-chemical model of disk atmospheres. 

Models of disk atmospheres are in an early stage of de- 
velopment. In part, this is because disks are complex 
systems, and many processes potentially play a role in 
determining the thermal-chemical properties of their at- 
mospheres. The relevant processes include those familiar 
from the study of molecular clouds and photodissocia- 
tion regions (e.g., UV and X-ray irradiation and their 
roles in heating, photochemistry, and ionization; neutral 
and ion-molecule reactions; and the freeze out and des- 
orption of molecules on grains). Other potentially im- 
portant processes are more unique to circumstellar disks 
and are generally less well understood, both in terms 
of how they operate and their impact on the thermal- 
chemical properties of disks. These include grain growth 
and settling, mechanical heating related to the magne- 
torotational instability, and radial and/or vertical trans- 
port of gas and solids. Models are developed to gain 
insight into how these processes affect disk atmospheres, 
with earlier studies having concentrated on different sets 
of processes and the interplay between them. 

In our previous work, we constructed a model of inner 
disk atmospheres that includes the thermal decoupling 
of the gas and dust (e.g., Glassgold et al. 2001, 2004; see 
also Kamp & Dullemond 2004; Nomura & Millar 2005), 



was based on X-ray heating and ionization, and included 
effects such as grain growth and accretion-related me- 
chanical heating (Glassgold et al. 2004). We used this 
model to study the gas-phase chemistry of water, includ- 
ing the effect of H2 formation on warm grain surfaces 
(Glassgold et al. 2009). Other studies have emphasized 
the effect of UV irradiation on the disk (e.g., Nomura & 
Millar 2005; Woitke et al. 2009; Woods & Willacy 2009; 
Bethell & Bergin 2009), included freeze-out and desorp- 
tion in the disk midplane or outer disk (e.g., Markwick 
et al. 2002; Woods & Willacy 2009; Walsh et al. 2010), 
and explored the effect of transport processes on disk 
chemistry (e.g., Ilgner et al. 2004, Semenov et al. 2006; 
Willacy et al. 2006; Heinzeller et al. 2011), processes 
that we have yet to include. In this paper, we build on 
our earlier work, implementing an improved X-ray ion- 
ization theory and an expanded chemical network and 
explore how variations in X-ray irradiation, accretion- 
related mechanical heating, and the 0/C ratio affect the 
thermal and chemical properties of the atmosphere. 

Our study complements other recent work on the 
chemistry of inner disks (e.g., Agiindez et al. 2008; 
Woitke et al. 2009; WiUacy & Woods 2009; Nomura et 
al. 2009; Kress et al. 2010; Markwick et al. 2002). The 
work of Agiindez et al. (2008) is perhaps the closest in 
spirit to our own study. In their study, they addressed 
the steady-state properties of the disk atmosphere, as- 
suming that the disk is strongly irradiated by FUV pho- 
tons and adopting an extended neutral chemistry net- 
work that includes carbon species. Our study comple- 
ments the Agiindez et al. study in that we include X- 
ray (but not UV) irradiation of the disk. While a more 
complete model would include UV irradiation, here we 
explore how well a model that excludes UV can do in 
accounting for the existing observations. 

2. DESCRIPTION OF THE MODEL 

We study the hydrocarbon chemistry of protoplane- 
tary disks with the steady-state thermal-chemical model 
developed by Glassgold, Najita & Igea (2004; hence- 
forth GNI04). This model is based on the D'Alessio et 
al. (1999) disk atmosphere for a generic T Tauri star of 
mass 0.5M,, blackbody temperature T, = 4000 K, and 
accretion rate 10~*M©yr^^. The disk is irradiated by 
stellar X-rays (which affect the temperature, ionization, 
and chemical structure of the gas), but not by far ultravi- 
olet (FUV) radiation. GNI04 used the dust temperature 
obtained by D'Alessio et al. and calculated the gas tem- 
perature on the basis of heating by collisions with dust 
grains. X-rays, and accretion-related mechanical heat- 
ing, balanced by a variety of line cooling processes and 
gas-dust collisions. As a result, the gas temperature is 
elevated over the dust temperature over much of the at- 
mosphere. 

We adopt a steady state treatment of the chemical ki- 
netics rather than a time-dependent one. This approach 
is supported in part by the warm gas temperatures in- 
ferred for the Spitzer molecular emission (300-1000 K), 
which generally imply short chemical timescales com- 
pared to transport timescales. The chemistry in GNI04 
was limited to 25 species for the elements H, He, C and O, 
and it contained 115 reactions. In a subsequent study of 
the water chemistry of protoplanetary disks, Glassgold, 
Meijerink & Najita (2009; henceforth GMN09) corrected 
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and expanded the chemistry to ~ 100 species and ~ 400 
reactions. 

GMN09 emphasized the role of H2 formation on warm 
dust grains, fohowing the theory of Cazaux & Tielens 
(2002, 2004; Cazaux et al. 2005) who showed that H2 
can form with moderate efficiency (< 0.2) on warm dust 
with temperatures up to ~ 900 K. Other recent studies of 
disk atmospheres have also employed this process in their 
models (e.g., Heinzeller et al. 2011; Woitke et al. 2009; 
Willacy & Woods 2009; and Woods & Willacy 2009). 
Enhanced H2 formation on grains appears to play an 
important role in producing significant warm molecular 
columns in models such as ours that include grain growth 
(Woods & Willacy 2009). If grain growth is not included, 
enhanced H2 formation on grains can lead to the over- 
production of warm molecular emission (e.g., Heinzeller 
et al 2011). 

In this paper we build on our previous models by em- 
ploying an expanded chemical network that includes hy- 
drocarbon, nitrogen and sulfur chemistry. The current 
model includes an enhanced theory of X-ray ionization 
that treats individually the nine most abundant heavy 
atoms (Adamkovics, Glassgold, & Meijerink 2011; hence- 
forth AGMll) in addition to C, N and O. The chemical 
network now contains ~ 115 species and ^ 1200 reac- 
tions. Even so, the molecular complexity is limited to 
molecules (and ions) of 6 atoms or less, with hydrocar- 
bons containing no more than two carbon atoms, as in 
C2II4. 

The reaction list contains neutral and ion-molecule re- 
actions. X-ray interactions, dissociative and radiative re- 
combination, and radiative association. We evaluated 
the reactions in our model using the primary literature 
and the critical evaluations of Baulch et al. (2005) for 
neutral reactions and Anicich (1993) for ion-molecule re- 
actions. We have also consulted the NIST Chemical Ki- 
netics Database (http://kinetics.nist.gov) as well as the 
UMIST (Woodall et al. 2007) and KIDA (Wakelam, 
2009) collections. 

Irradiation by stellar and/or interstellar FUV is not in- 
cluded in the present model, an assumption that would 
apply when the disk is well shielded from UV (e.g., by 
an intervening funnel flow). This assumption may also 
be more generally valid in the context of heating. Ear- 
lier work (e.g., Nomura et al. 2007) has shown that UV 
heating due to the photoelectric effect is significantly re- 
duced relative to other heating processes when grains 
have grown and settled to the extent that has been de- 
duced for T Tauri disks (Furlan et al. 2006). This, and 
the uncertain PAH abundance of T Tauri disks, mitigates 
our neglect of UV heating, a topic that we return to in 
section 4. More seriously, UV photodissociation likely 
affects the molecular abundances as we also discuss in 
section 4. 

One obvious difference between our X-ray irradiated 
model and others that focus on UV irradiation is that in 
our model molecules are destroyed by H+ and Hc+. In 
the presence of UV irradiation, photodissociation would 
also play a significant role in destroying molecules. In 
a similar vein, UV irradiation will also alter the ion- 
ization and temperature structure from that considered 
here. Another difference between UV and X-ray irradi- 
ated models is that X-rays are also able to ionize species 
that are unavailable to UV photons. For example, al- 



though FUV photons are unable to directly ionize N and 
N2 , ion- neutral reactions of N2 with and He"*" (which 
are generated by X-rays) can activate the nitrogen chem- 
istry. 

We include the formation of H2 on dust grains in the 
same way as GMN09 (their eqs. 14 and 15), following the 
theory of Cazaux & Tielens. At the suggestion of Bai 
& Goodman (2009), Cazaux & Tielens (2010) have since 
corrected their calculations, but the changes are only of a 
quantitative nature. We ignore the possibility of forming 
other molecules on grains. GMN09 found that ignoring 
H2O formation on grains did not much affect the distri- 
bution of gaseous H2O. We also ignore cooling by H2O 
rotational transitions as well as heating due to UV ab- 
sorption by H2O. Water line emission has been measured 
to be a significant coolant for disk atmospheres (Pontop- 
pidan et al. 2010). These limitations will be remedied in 
future reports. 

As discussed in more detail in AGMll, the abundances 
of the heavy elements play an important role in the 
chemistry of the disk atmosphere, because they are the 
dominant charge carriers below the atomic to molecu- 
lar transition. Measurements are not available for the 
absolute gas-phase abundances of the elements in thick 
molecular clouds, not to mention protoplanetary disks. 
The only abundance determinations come from UV ab- 
sorption line spectroscopy of relatively diffuse interstellar 
clouds {Ay ^ 2 mag) obtained from observations made 
with satellite observatories. The measured abundances of 
the heavy elements vary from cloud to cloud, and we con- 
tinue to use those abundances reported for the very well 
studied cloud towards ( Oph (Savage & Sembach 1996); 
as shown in Table 1. The last column shows the deple- 
tion factors, defined as the ratio of solar to interstellar 
abundances, with the former coming from Table 1 of the 
review of solar abundances by Asplund et al. (2009). As 
discussed in this paper and in Jenkin's (2009) review of 
interstellar abundances, both sets of values have signifi- 
cant uncertainties, in some cases factors of a few. Jenkins 
highlights the case of S, where the UV observations sug- 
gest essentially no depletion in interstellar clouds, despite 
the fact that this element is seen in meteorites. 

Although we use the depleted interstellar abundances 
in Table 1 for the preliminary modeling purposes of this 
paper, we recognize that the abundances of the heavy el- 
ements in diffuse interstellar clouds may well have been 
altered in dense clouds and by the formation and evolu- 
tion of protoplanetary disks. The determination of rela- 
tive abundances in disks remains a primary goal of ob- 
servations of the emission line fluxes of their atoms and 
molecules. As discussed in Section 1, we will examine 
one particular modiflcation of the elemental abundances, 
the variation of the oxygen-to-carbon ratio (0/C) that 
may be induced by the adsorption of water on solid bod- 
ies and the transport of those bodies between large and 
small radii. 

The molecular abundances in the inner disk atmo- 
sphere are affected by the level of X-ray irradiation, the 
strength of the accretion-related mechanical heating, and 
the oxygen-to-carbon ratio (0/C). Of course they also 
change with the radial distance from the star. We have 
investigated these effects by carrying out variational cal- 
culations about a "reference case" defined in Table 2. 
The X-ray variations are parameterized by the X-ray 
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luminosity Lx and the X-ray temperature Tx, and the 
mechanical heating by the parameter ah in Eq. (12) in 
GNI04, 

9 

Tacc = -a;ipc^ri, (1) 

where p is the local mass density, c is the isothermal 
sound speed, ft is the angular rotation speed and ah 
is a phenomenological parameter specifying the strength 
of the heating. As discussed by GMN09 in connection 
with their Equation 1, the grain surface area per unit 
volume is proportional to the dust-to-gas ratio of the 
disk, pdust/p, and the inverse of the geometric mean (cg) 
of the minimum and the maximum grain radii in an MRN 
dust size distribution (n(a) oc a~^-^). 

Table 2 gives the range of parameter variations we con- 
sidered. As in our previous papers, the dust-to-gas ratio 
was fixed at 0.01, and we considered values of Ug that 
are 32 and 320 times larger than the usual interstellar 
value of 0.0224 /xm, i.e., 0.707 and 7.07 /xm. This range 
is consistent with the dust depletion factors (^100 1000) 
that arc found for the upper disk layers of T Tauri stars 
based on their mid-infrared colors (Fiulan et al. 2006). 
We adopted Lx and Tx values that are typical of T Tauri 
stars {Lx = 10'^^ - 10^^ ergs"^; Feigelson et al. 2005) as 
well as lower values of Lx down to lO^^ergs"^ in order 
to illustrate the role of the X-rays. 

The parameter ah was varied so that accretion heat- 
ing was either non-dominant (ah = 0.01) or dominant 
{ah = 1) as a heat source. The lower value is motivated 
by the viscosity parameter for T Tauri disks, which has 
been estimated as a ~ 0.01 as a vertical average in a 
steady state disk (e.g., Hartmann et al. 1998). The larger 
value is motivated by 3D MHD calculations that suggest 
that turbulent accretion disks driven by the magnetoro- 
tational instability experience greater energy dissipation 
at the disk surface than in a classical "a-disk". A re- 
cent calculation by Hirosc & Turner (2011) suggests that 
the energy dissipation rate rises toward the disk surface, 
reaching values that correspond to ~ 1. 

The range of radii we considered covers the inner disk 
region probed by Spitzer (the inner few AU) and the 
larger radii relevant to Herschel (out to 40 AU). We con- 
sidered a range in the 0/C ratio in order to explore 
the dehydrated to supcr-hydrated regimes described in 
Ciesla & Cuzzi (2006). Those models indicate that the 
(vertically averaged) water vapor concentration can in- 
crease by up to an order of magnitude (e.g., their Fig- 
ure 1), or decrease by more than an order of magnitude 
(e.g., their Figure 3) in the inner few AU region. We 
translated this into an enhancement or decrement in the 
oxygen abundance by up to an order of magnitude rela- 
tive to the standard case. We describe the results of our 
model calculations for the reference case in Section 3.1 
and the results for the variations in Section 3.2. 

3. RESULTS 

3.1. Reference Case 

Figure 1 shows the vertical variation of some of the ba- 
sic physical quantities of the disk atmosphere for a range 
of radial distances in the inner disk between 0.25 AU and 
20 AU. In all three panels, the abscissa is the vertical 
column density A^^h, integrated downward into the disk 
from above. The top panel shows the ionization rate mul- 



tiplied by the inverse-square dilution factor. The similar- 
ity of the profiles at the different radii shown illustrates 
the usefulness of characterizing disk properties as a func- 
tion of vertical column. Note that Compton scattering, 
which is not treated here, begins to become important 
for column densities approaching the maximum value in 
the figure A^h = lO^'^cm"^. 

The temperature variation with column density (Fig. 1; 
middle panel) shows a particiilarly rapid drop at verti- 
cal column densities Na ^ 0.6 — 1.5 x lO^^cm"^. The 
temperature is high, in the range 4000-5000 K, at col- 
umn densities below this value, but it then drops into 
the warm 200-2000 K range at larger column densities. 
The atmosphere remains warm until collisions with dust 
grains bring the gas temperature down to the dust tem- 
perature; this occurs somewhere in the column density 
range A^h = 10^^ — lO^^cm"^, with the value depend- 
ing on radius. The 200-2000 K region (henceforth the 
"warm atmosphere") is a region of high density with 
Tin = 3 X 10^ — 3 X 10^° cm~^ in the reference case. 
As shown in the bottom panel of Fig. 1, the ionization 
fraction decreases by factors of 100 where the gas tem- 
perature drops abruptly and the transition from H to H2 
occurs. 

Many molecular abundances also increase by large 
amounts as the temperature and ionization fraction drop, 
led by the atomic to molecular hydrogen transition shown 
in the bottom panel of Figure 1. The change from a 
primarily atomic to a primarily molecular medium, and 
from higher to lower levels of temperature and ioniza- 
tion, are the manifestations of a thermal- chemical tran- 
sition. This transition is mediated thermally by the 
onset of strong molecular line cooling that causes the 
sharp decrease in temperature. The chemical transi- 
tion to molecules results from the decrease in the ion- 
ization level of and He^, which destroy molecules, 
and the suppression of chemical temperature-dependent 
reactions that destroy them. 

Thermal-chemical transitions are commonly found 
among models in the literature, although the details of 
the transition may differ depending on the irradiation 
field and other model properties. The sharpness of the 
transition is one property that can vary depending on 
the model parameters, and sharp transitions like that of 
the (X-ray irradiated) reference model may appear in UV 
irradiated models as well. The UV irradiated model of 
Woitke et al. (2009) is a good example. The electron frac- 
tion in that model varies sharply, from 10"** to 10~^ at 
z/r = 0.15 over a radial range from ~ 0.5 to 10 AU. The 
specific details differ from our model, although the possi- 
bility of sharp transitions is a common clement. Models 
that include both UV and X-rays also find results similar 
to ours, e.g., a recent calculation by Woitke (2011, per- 
sonal communication) has a similar temperature struc- 
ture. Nomura et al. (2007) also find a sharp temperature 
decline when X-rays and UV are considered and grains 
have grown to large size. There is also diversity among X- 
ray irradiated models: we find sharper transitions when 
mechanical heating is included, and shallower transitions 
when mechanical heating is weak (Glassgold et al. 2004, 
2009; see also the following section). 

The top two panels of Figure 2 show the transitions 
of carbon, oxygen and nitrogen to molecules for the ref- 
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erence model at radii 0.25, 1.0, and 2.0 AU (Figure 2a), 
and 4.0, 10, and 20 AU (Figure 2b). The abundances 
in this figure are defined relative to that of all hydrogen 
nuclei, and the region below the main thermal-chemical 
transition near TVh ~ 1 — 3 x 10^^ cm~^ is emphasized. 
Focusing for the moment on the results for r = 1 AU, 
we find that, like the transition from atomic H to H2, 
the transition from atomic carbon to CO is complete, 
and CO remains abundant down to the largest vertical 
columns. In contrast, the conversion of all of the residual 
oxygen (that is not in CO) into H2O occurs in the warm 
atmosphere (TVr = 1 - 3 x lO^^cm'^), but the H2O 
abundance declines at larger columns where the gas is 
no longer warm, with O2 and SO2 as well as atomic O 
becoming the major reservoirs of the residual oxygen. 

Several other recent thermal-chemical calculations also 
study H2O synthesis in disk atmospheres. In the model 
of Woitke et al. (2009), water is present in the disk at- 
mosphere at 1 AU although at much lower column densi- 
ties (< 10^^ cm~^ above 200 K) than found here. In the 
calculation of Willacy & Woods (2009), a much larger 
column of warm water (^^ 10^® cm~^) is present and at 
higher temperature (> 1000 K). 

As shown in the second panel (from the top) of Figure 
2, the transformation from atomic to molecular nitro- 
gen becomes complete in the warm zone, and atomic N 
remains fairly abundant throughout much of the inner 
disk. The bottom panel shows the molecular ions gener- 
ated by X-ray ionization. Although they are important 
for the chemistry, they do not dominate the ionization 
balance, which is accomplished in this model by heavy 
atomic ions such as Si^, Fe"*" and Mg"*". However, this 
conclusion is sensitive to the assumed elemental abun- 
dances, as discussed by AGMll. 

We find that a thermal-chemical transition from hot 
and atomic to warm and molecular always occurs in X- 
ray irradiated disks, at least out to 40 AU, although the 
details of the transition (such its location, sharpness, and 
other properties) are sensitive to the input parameters of 
the chemical and thermal model. We discuss this topic 
further in Section 3.2, where we describe variations about 
the reference model. In this section we focus on the 
molecular abundances found in the reference case. 

The physical and chemical properties of the disk atmo- 
sphere that are reported here are consistent with those 
from our previous models (e.g., GNI04, GMN09). The 
present work does have some quantitative differences 
with our earlier work. These are associated with the 
many improvements in the physical-chemical model and 
the computational program employed here. 

It is no accident that the transition to from atomic 
C to CO follows closely the temperature and ionization 
transitions in Figure 1, since CO is a dominant coolant 
in the thermal-chemical transition. The attainment of 
maximal CO abundance occurs relatively high in the disk 
atmosphere through the action of warm neutral radical 
reactions, as discussed in GNI04 (Sec. 3) and GMN09 
(Sec. 2.2). The formation of the OH radical is a key step 
in the synthesis of CO and H2, through the endothermic 
reaction, 

O + H2 ^ H + OH; (2) 
the rate coefficient recommended by Baulch et al. (1992) 



k = 8.5 X 10-2" T^-^^ e-3i63/Tcm3g-i^ 

Water is then formed by an exothermic reaction with a 
modest barrier, 

OH + H2^H20 + H; (4) 
its rate coefficient is 

fc = 8.5 X 1.70 X 10-16 Ti f^° e-i660/Tcm3g-i_ (5) 

By contrast, CO is formed from OH by a fast neutral re- 
action without a significant barrier. Despite the barriers 
in Eqs.[3]and[5l water and CO are formed efficiently by 
the above reactions at temperatures as low as ~ 250 K. 

This situation is not duplicated in the case of hydride 
formation for the two other abundant and chemically in- 
teresting heavy atoms, C and N. In reactions analogous 
to Eq. 131 with O replaced by C and N, the thermal bar- 
riers are now ^ 12000 K instead of ^ 3000 K. These bar- 
riers suppress the formation of CH and NH by reactions 
with H2, and the inverse reactions, in which atomic H de- 
stroys CH and NH, are both fast. In our X-ray irradiated 
models, the abundance of H decreases slowly below the 
thermal-chemical transition, so both suppressed forma- 
tion (via reactions with H2) and rapid destruction (by H; 
i.e., reactions analogous to those in Eqs. 1 and 3) imply 
that neutral production of carbon and nitrogen hydrides 
along the lines of oxygen-hydride synthesis is inefficient. 

Nonetheless, the warm neutral chemistry that is re- 
sponsible for most carbon and oxygen going into CO and 
H2O plays an important role in the synthesis of other 
molecules. This is shown in Figure 3, which is a highly 
schematic chemical network that shows how fast neu- 
tral reactions of OH with N, O and CO generate other 
molecules such as CO2, N2, and SO2 (not shown) and 
more complex molecules that will be discussed below. 
Thus the OH radical is the centerpiece of the warm chem- 
istry of protoplanetary disk atmospheres. 

As shown in Figure 3, the production of NO by the 
fast reaction, 

OH + N^NO + H fc ~ 5 X 10-"cm3s-\ (6) 

leads to the formation of N2 

N0 + N^N2 + fc -4 X lO-^cm^s-^ (7) 

The nitrogen hydrides can be synthesized by successive 
hydrogenation of the N+ ion, produced by X-ray ion- 
ization, via successive reactions with H2, followed by 
dissociative recombination of the hydride ions. lon- 
hydrogenation becomes weak beyond NHs"'', but proton 
transfer from molecular ions to neutral hydrides will pro- 
duce hydride ions up to NH4'. Thus a combination of 
ion-molecule and neutral reactions initiated by X-ray ion- 
ization can synthesize ammonia. A related sequence of 
neutral reactions starting with NO leads to the synthesis 
of HCN, 

NO + C^CN + fc - 3 X 10-"cm3s-\ (8) 
CN+H2 ^ HCN+H fc(300K) - 2.0xl0-i'^cm3s-\ (9) 

Although many reactions have been omitted from this 
brief description of the nitrogen chemistry for purposes 
of clarity, the middle panel of Fig. 2 shows that the main 
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C CO and O — > H2O transitions arc accompanied 
by the almost full conversion of N to N2. The transfor- 
mation from N to N2 is mediated by neutral reactions 
as indicated in Fig. 3. It is accompanied by significant 
warm columns of other nitrogen-bearing molecules, such 
as NO, NH3, HCN and N2H+. 

The integrated treatment of both neutral and ion- 
molecule reactions is a feature of the underlying chem- 
ical model, with sonic; molecular species demonstrat- 
ing a greater sensitivity to X-ray irradiation than other 
species. The X-rays play a critical role in the synthesis of 
NH3 and N2H+. Since FUV radiation is unable to ion- 
ize N and N2 (which have ionization potentials of 14.5 
and 15.6 eV, respectively), the X-rays are important in 
activating the nitrogen chemistry. 

The hydrocarbon chemistry proceeds similarly to the 
nitrogen chemistry. The number of hydrogen atoms in 
a hydrocarbon molecule or ion can be increased by re- 
actions with H2 or decreased by reactions with H. Neu- 
tral hydrocarbons can be ionized by proton transfer from 
abundant molecular ions, such as Hj, HCO+ and H3O+, 
and neutralized by dissociative recombination with elec- 
trons; their charge is also affected by charge transfer pro- 
cesses. These types of reactions are believed to be oper- 
ative for the hydrocarbon families of interest, those with 
one or two (;a.rbon atoms. A key point is that double- 
carbon molecules are synthesized from those with just 
one carbon atom by insertion reactions with C and C+, 
as illustrated by the simplest examples, 

C + CH„ ^ C2H„_i + H (n = l-4), (10) 

C+ + CH„ ^ C2H+_i + H (n = l-4), (11) 

or by molecular ions as suggested by Woodall et al. 

(2007) 

cut + C ^ C2H+ + H2 and CH^ + C-^ C2H+ + H. 

(12) 

These highly simplified remarks about hydrocarbon 
chemistry are closely aligned with the general discussion 
of Agiindcz, Ccrnicharo & Goicocchea (2008). In their 
model variation that is closest to our model, these au- 
thors use a photodissociation code to calculate the chem- 
ical structure of the surface region of the D'Alessio et 
al. (1999) disk model. They assume that the gas and 
dust temperatures are the same and that the dust is in- 
terstellar, and they consider stellar as well as interstel- 
lar UV radiation but ignore X-rays. Despite the many 
differences with their model, we draw the same general 
conclusion: hydrocarbons are predicted to be present in 
the top layers of protoplanetary disk atmospheres. There 
are of course quantitative differences between the results 
of Agiindcz et al. and those reported here. We empha- 
size warm molecular column densities that arise in re- 
gions with T > 200 — 300 K, because they are related 
to measured fluxes. Agiindez et al. do not calculate gas 
temperatures, and they report total vertical column den- 
sities in the disk surface layer, almost all of which is cool 
gas (< 300 K) at lAU. 

Willacy & Woods (2009) also synthesize a signficant 
column of warm HCN. At 1 AU they find a column of 
^ 4 X 10^^ cm~^ at a temperature > 1000 K. As they 



note, the much lower columns of C2II2 that result in their 
models ('-^ 9 x 10® cm^^) are likely due to the omission 
of reactions involving H2 with C2 or C2H. 

Despite the fact that CO tends to take up most of 
the available carbon, the hydrocarbon chemistry is still 
active after the main transition. This is a characteristic 
feature of the X-ray chemistry, where X-ray generated 
He+ ions produce C+ ions when they destroy CO by the 
ion-molecule reaction, 

He+ CO C+ O He. (13) 

Similar reactions, in which He"*" generates atomic ions 
from other molecules, are in fact quite common, notably 
for N2 and 002- In the latter case 0+ is also produced. 

In general, we flnd that both neutral and ion-molecule 
reactions contribute to the production and destruction of 
molecular species, with the relative importance of the two 
kinds of reactions varying with position in the disk as well 
as with species. Neutral reactions are more important in 
the warm disk surface, while ion- molecule reactions play 
a significant role in the cooler regions below. The species 
shown in Fig. 2 offer many challenges for the detection of 
molecules in the inner parts of protoplanetary disks, and 
some have already been observed, as discussed in Sec. 1. 

Wc briefly summarize the results for the reference 
model at 1 AU, describing our results for the more 
abundant molecules. To do this, we take advantage of 
the monotonically decreasing run of temperature with 
increasing vertical column density (Fig. 1) and inte- 
grate the species' volumetric densities downward from 
the maximum (2000 K) to minimum (200 K) tempera- 
ture of the warm atmosphere. This corresponds to the 
region just below the main thermal-chemical transition 
discussed in Section 2 and roughly to the region that 
can produce observable line emission in disks seen partly 
face-on. We often refer henceforth to these integrated 
columns simply as "warm columns" . 

• H2O forms in the same region as CO and reaches a 
comparable abundance, ~ 10~*. The water abundance 
declines below 4 x 10^"'^ cm~^ before recovering its high 
abundance approaching lO^^'cm^^. The resulting col- 
umn density of warm water is ^ 5 x 10^^ cm^^. In the 
region of reduced water abundances, several other species 
take up the available oxygen. Together with water they 
(i.e., O2 and SO2) would enhance the cooling above that 
provided by O and CO, the two species that are con- 
tribute to the cooling in the present calculation. 

• HCN is less abundant, reaching a maximum abundance 
of ~ lO^'^ in the warm atmosphere. The total column of 
warm (> 200 K) HCN is ~ lO^^cm'^. 

• C2H2 is less abundant than HCN, reaching a maximum 
abundance of 10~® in the warm atmosphere. The column 
of warm C2H2 is ^ 10^'^ cm^^. 

• CO2 has an abimdance in the range 10"^ — 10~® 
in the warm atmosphere; its warm column density is 
~ 5 X 10i''^cm-2. 

• HCO is more abundant than CO2, with a warm column 
density of ~ 1.4 x 10^^ cm~^. The warm column density 

of HCO+ is - 1.4 X 10^3 cm-2 

• NH3 reaches its maximum abundance of ~ 10^^ just 
below the warm region; its warm column is ^ 3 x 
10^5 cm-2 

While H2O and some other molecules such as HCN and 
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C2H2 form high up in the atmosphere and then dechne 
in abundance, some species increase further down (e.g., 
CH4), or at least do not decrease (NH3; Fig. 2). CH4 
has a maximum abundance in the middle of the warm 
region where T ^ 800 K. Just below the warm region, it 
begins to recover from its abundance low (< 10"^") and 
becomes abundant at large depths: a;(CH4) ^ 10^'' for 
A^H ^ 10^'^ cm~^. We find that the chemical as well as the 
temperature and ionization structure of our disk model 
is highly stratified. Such a stratified structure is unlikely 
to be altered by turbulent mixing unless the transport 
timescale is unusually short (e.g., Heinzeller et al. 2011). 

3.2. Variations on the Reference Case 

In the previous section we presented the abundances 
of major species obtained for the reference case defined 
in Table 2. Even for fixed stellar and disk properties, the 
model is characterized by several parameters: the me- 
chanical (accretion) heating, the dust surface area, the 
C/0 ratio and the radial distance. Here we investigate 
how the variation of these parameters affects the predic- 
tions of the model. 

One important uncertainty is the efficiency of H2 for- 
mation on grains, as discussed at the beginning of Sec- 
tion 2. Without it, i.e., with only gas phase formation of 
H2, the abundances of molecules other than CO (H2O, 
CO2, HCN, C2H2) are much reduced in the warm atmo- 
sphere, with the result that their warm columns are very 
smah (see also Agiindez et al. 2008; GNI04). As noted 
in Section 2, theoretical calculations do support H2 for- 
mation on warm grains, although there is still consider- 
able uncertainty about the efficiency of this process (e.g., 
Cuppen et al. 2010). The formation rate also depends on 
the surface area of the grains, which we parameterize in 
terms of Og (defined in Sec 2). The grain surface area in 
the reference model is smaller than the ISM value by a 
factor of ^ 30. 

Figure 4 compares the thermal-chemical transition of 
the reference case (black curves with ag — 0.707/xm, 
ah = 1.0) and two other cases in which Og is increased 
by a factor of 10 (red curves) and au is reduced by 100 
(blue curves). The red curves correspond to a reduction 
in the grain surface area per hydrogen nucleus by a fac- 
tor of ^ 300 relative to interstellar grains, and the blue 
curves have essentially no mechanical heating. In the 
former case, the transition is shifted significantly deeper 
down into the disk, from 1 x lO^^cm"^ to 5 x lO^^cm"^, 
with a corresponding increase in volumetric density. Be- 
cause the dust cooling rate is also reduced by a factor 
~ 300, the depth of the warm layer is extended, and 
the warm molecular columns, given in Table 3, are all 
increased, with the exception of OH. When mechani- 
cal heating is essentially eliminated, the transition is 
smoother and shifted upwards in the atmosphere to the 
region 2 x 10^°cm~^ — 5 x 10^°cm~^. The associ- 
ated warm region is also shifted upwards, and the warm 
molecular columns are decreased by a factor of five to 
ten, with the exception of OH, which is increased by a 
factor of about three. 

From this discussion of we can conclude that warm 
molecular columns increase with increasing dust size (i.e., 
reduced dust surface area; parameter Og) and increasing 
mechanical heating (parameter au)- Both effects are ac- 



companied by a downward shift in the location of the 
atomic to molecular transition. The warm region starts 
at the transition, and larger warm columns (which favor 
detection) are generated as the transition moves down- 
ward into regions of higher density. Table 3 gives the 
warm columns for some of the more abundant species 
that have either have already been observed in T Tauri 
disks or might be candidates for future detection. 

X-rays are the only source of external radiation in 
these model calculations. While varying Tx from 1 keV 
to 5 keV did not much alter our results, variations in 
Lx have a much larger effect. Table 4 shows how the 
warm column of selected molecules and ions depends on 
Lx at lAU for a large range of values of Lx when the 
other model parameters are kept the same as in the ref- 
erence case. From Chandra observations of young clus- 
ters, Preibisch & Feigelson (2005) found that the X- 
ray luminosity function depends mildly on cluster age 
and, within clusters, on the stellar mass range. For 
the Orion Nebular Cluster, the median K star X-ray 
luminosity is ix = 3 x 10^" erg s"^, with 90% of the 
stars having Lx > 4 x 10^^ erg s""'^ and 10% having 
Lx — 1-3 X lO^^ergs"^. In Table 4 we can see that the 
warm columns of ions are especially sensitive to the X-ray 
luminosity and, as might be expected, increase dramat- 
ically as Lx increases. Even for X-ray luminosities as 
low as Lx = lO^^ergs"^ (i.e., several orders of magni- 
tude lower than the X-ray luminosities of typical T Tauri 
stars), the X-rays are a significant source of ionization, 
and could influence the chemistry of T Tauri disks by 
producing a variety of molecular ions by proton transfer 
and by destroying neutral molecules by reactions with 
H+ and He+. At these low X-ray ionization rates, the 
stellar UV radiation field (not treated here) would also 
be particularly important. 

But simple neutral molecules change little as Lx goes 
from 10^^ to 10^^ ergs~^. The reason is that mechanical 
heating associated with the MRI theory of viscosity dom- 
inates the heating of the disk atmosphere, and molecular 
synthesis in the warm regions is built on temperature- 
sensitive neutral reactions, as discussed in the previ- 
ous section. In the presence of strong accretion heat- 
ing, the X-rays also have little effect on the location of 
the main thermal-chemical transition defined by CO and 
H2. For example at i? = 1 AU, the transition shifts from 
6.5 X lO^"^ cm-2 for Lx = 10^^ erg s'^ to 1.0 x lO^^ cm'^ 
for Lx — lO^^ergs"^, while for Lx < lO^^ergs"^ there 
is no change at all. Nonetheless, molecular ions still play 
a role, sometimes in forming neutrals by dissociative re- 
combination and sometimes in destroying them by pro- 
ton transfer or charge transfer reactions. 

By way of contrast, radicals are sensitive to X-rays, as 
exemplified by OH and NO, both of which increase by 
about 3,000 when Lx is increased from lO^^ergs"^ to 
j^Q3igj.gg-i jj^ Qyj. model, the steady state abundance 
of OH in warm regions is mainly determined by its pro- 
duction by recombination of H3O+ and destruction by 
neutral reactions. NO follows OH, since it is produced 
from OH via Eq. [6]and destroyed by atomic N via Eq. [71 
A similar explanation applies to SO and SO2, which are 
moderately sensitive to X-rays. It is of considerable in- 
terest that the warm columns of the hydrocarbons CH4 
and C2H2 decrease with increasing Lx. This trend is 
stronger for CH4, very likely due to the fact that it is 
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produced by dissociative recombination of the molecular 
ion CH^ which is produced by the radiative association 
of CH+ an H2. 

Variations in the 0/C ratio can have a dramatic effect 
on the abundances in the warm disk atmosphere. Ta- 
ble 5 gives examples for the reference case at 1 AU. The 
oxygen abundance has been varied at fixed carbon abun- 
dance so that the 0/C ratio varies from 0.25 to 5.0; the 
reference case value is 2.5. Below 0/C ~ 1, the reduction 
in the available oxygen leads to a significant decrease in 
the warm column of CO and also of CO2. Pure oxygen 
species such as O, OH, H2O, O2, HsO"*" show very large 
decrements, especially below 0/C = 1, as do closely re- 
lated oxides such as SO2 and NO. At the same time, pure 
carbon species like C2H2 and CH4 increase drastically. 
Nitrogen species also change, but by large amounts only 
when they contain carbon (e.g., HCN) or oxygen (e.g., 
NO). One noteworthy result of this analysis is that the 
warm column of Ne"*" increases rapidly with decreasing 
0/C ratio because it is efficiently destroyed by H2O. 

Ratios of warm columns derived from molecular line 
flux measurements are often used in making preliminary 
interpretations of the observations. Our results show 
that variations in the 0/C ratio can cause significant 
changes in the warm column ratios. Because H2O de- 
creases more rapidly than CO with decreasing 0/C, the 
H2O/CO ratio becomes very small for small 0/C < 1. 
By way of contrast, the HCN/H2O ratio is a sensitive 
decreasing function of the 0/C ratio and becomes large 
as 0/C approaches unity. The same is true for the 
C2H2/H2O ratio, but even more so because the warm 
column of C2H2 is a much more sensitive function of 
0/C than HCN. Thus these ratios can be considered as 
signatures of enhanced carbon chemistry. 

Looking at how the overall abundance patterns vary 
with disk radius, we find that the abundance patterns at 
1 AU are similar at smaller and larger radii, as demon- 
strated in Figure 2 for oxygen molecules (top panel), ni- 
trogen molecules (middle panel) and molecular ions (bot- 
tom panel). Some species show large variations in their 
vertical column density distributions as a function of ra- 
dius, e.g., CO2. Our results for 40 AU are not plotted 
because they are qualitatively similar to those for 20 AU. 
Figure 5 and Table 6 show the radial variation of warm 
columns (i.e., the regions of the disk where line emis- 
sion can be efficiently generated) for the reference case. 
The radial range in Figure 5 and Table 6 extends only to 
10 AU because the disk atmosphere at 20-40 AU does not 
have significant column density in the 200-2000K range 
(Fig. 1). Note that the gas column densities shown in 
Figures 2 and 5 and Table 6 do not include the effect of 
freeze-out, which would have a greater impact in the disk 
midplane. 

We see that the warm columns are large at very small 
radii, decrease rapidly inside 1 AU, and then decrease 
more slowly beyond. This behavior is typified by H2 itself 
and is shared by many of the more abundant molecules, 
e.g., CO, H2O, O2, NH3, and HCN. The water column 
increases significantly at radii < 0.5 AU because the neu- 
tral reactions that produce water are effective throughout 
the warm atmosphere (because of the high gas tempera- 
ture and density in this region) and lead to a high abun- 
dance of water throughout (Fig. 2). In contrast, at 1 AU, 



the water abundance drops at intermediate heights (be- 
low iVn = 5 X 10^^ cm^^) before recovering deeper down, 
leading to a lower warm column than at 0.25 AU (Fig. 5). 

It is interesting that, in the top panel of Figure 5, the 
large warm column of atomic H changes by no more than 
a factor of two from 0.25 AU to 10 AU. The reason for 
this is that the basic OH chemistry operative in the warm 
region (Eqs.[2]and|4]) releases atomic H as atomic oxygen 
is processed into water. The abundance of atomic H 
is important for the synthesis of hydrocarbon molecules 
because many of the reactions that hydrogenate carbon 
are reversed by atomic H reactions. 

Many of the heavier molecules show similar patterns 
because they share a common (warm neutral radical) 
chemistry and because the physical properties of the 
warm region change slowly with disk radius. For ex- 
ample, the ionization parameter has the value C/'^h ~ 
3 X 10~^''cm^s~^ to within a factor of two over radii from 
0.25 to 20 AU (a factor of almost 100) because the inverse 
square dilution of the X-ray ionization rate (top panel of 
Fig. 1) is largely compensated for by the decrease of the 
density with radius. Consequently the electron fraction 
at the top of the warm region is almost independent of 
radius, with a value close to = 2 x 10"^. This approx- 
imate invariance applies to the ionization of the warm 
region, but not to the molecular abundances which are 
affected by density and temperature (and processes such 
as freeze out) as well as ionization. The temperature 
does remain close to 1000 K at the top of the warm col- 
umn for small radii (middle panel of Fig. 1), but it then 
decreases significantly beyond 4 AU. Of course the den- 
sity at a given vertical column density decreases by more 
than 1,000 as the radius increases from 0.25 to 20 AU. 
These changes also cause the thermal chemical transition 
to move slightly upward in the atmosphere with increas- 
ing radius, and this also tends to lead to a reduction in 
warm column densities. The decreasing density of the 
warm atmosphere as a function of radius may play a role 
in limiting the radial extent of emission from molecular 
transitions with high critical densities. 

4. DISCUSSION 

4.1. Inner Planet Formation Region 

Our model results (Section 3) have potential implica- 
tions for the understanding of the inner and outer planet 
formation regions of the disk, i.e., the regions roughly 
within and beyond the "snow line" where water vapor 
would condense on grains. The inner planet formation 
region, discussed in this section, is probed by the molec- 
ular observations made with Spitzer and ground-based 
near-infrared spectroscopy. The outer planet formation 
region is discussed in the next section. 

4.1.1. Warm Molecular Columns 

As we found in GMN09, the warm atmosphere in our 
models achieves temperatures comparable to that of the 
Spitzer molecular emission from T Tauri disks (~300- 
1000 K; Carr & Najita 2008, 2011; Salyk et al. 2008, 
2011). These temperatures are maintained over signif- 
icant column densities at disk radii out to ~ 4AU (Fig- 
ure 1, middle panel). Thus the temperature structure of 
our disk atmosphere is consistent with the interpretation 
that the Spitzer molecular emission arises from the inner 
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few AU of the disk. 

Some molecular species appear quite sensitive to the 
parameter variations we considered, with their warm 
molecular columns varying significantly over the param- 
eter range we explored. As one example, the variation 
in grain size and mechanical heating alone can produce 
warm column variations of ^ 30 in H2O and HCN (Table 
3). This sensitivity may help to account for the factor 
of ~ 30 range in the H2O and HCN emission strengths 
among T Tauri disks (Pontoppidan et al. 2010; Salyk 
et al. 2011; Carr & Najita 2011). The role of grain set- 
tling in producing enhanced warm columns may also help 
to account for the finding by Salyk et al. (2011) of the 
enhanced detectability of molecular emission from more 
settled disks. 

The warm molecular columns we find for some species 
are also quantitatively similar to those inferred from the 
Spitzer data. The line-of-sight column density of the wa- 
ter emission detected with Spitzer (~ lO-'^^cm"^), which 
has an associated projected emitting area of Trr^ where 
r ~ lAU (Carr & Najita 2011; see also Salyk et al. 
2011). is similar to our (vertical) warm column of wa- 
ter at lAU (10^^ - 5 X lO^^cm-^; Tables 3, 4, and 5). 
The much smaller warm columns achieved at larger disk 
radii (< 10^^ cm^^ at 10 AU in the reference case) help 
to account for the inference of AU-sized water emission 
regions from the Spitzer spectra, although perhaps not 
completely (see Section 4.2). 

Larger warm columns are achieved at small disk radii 
(~ lO^^cm"^ at 0.25 AU in the reference case; Table 
6 and Figure 5). These are reminiscent of the water 
columns reported at such radii based on K- and L-band 
spectroscopy of disks (~ 10^° — 10^^ cm~^; Carr et al. 
2004; Najita et al. 2009; Doppmann et al. 2011). How- 
ever, the temperature of the observed emission is warmer 
(~ 1500 K) than most of the warm column (2000 200 K) 
in the model. Thus, while the empirical situation is that 
strong near-infrared water emission is rare (e.g., Najita et 
al. 2007), this does not necessarily mean that water is un- 
der abundant at small disk radii in most disks. Our mod- 
els suggest that large columns of water are commonly 
present at small disk radii, but the water is typically not 
warm enough to produce strong water emission in the 
near-infrared. 

We can also compare our results to those derived for 

the 4.7 /im CO fundamental emission from T Tauri stars, 
which has a characteristic temperature of ^500 1000 K, 
line-of-sight column density ^ 10^* cm^^, and projected 
emitting area of radius ~ 0.15 AU (Salyk et al. 2011; 
Najita et al. 2003). Our reference model at 0.25 AU pro- 
duces a vertical CO column of 5 x 10^^ cm^^ over the ver- 
tical temperature range 2000 500 K, in reasonable agree- 
ment with the observed values. The decrease with radius 
of both the density of the atmosphere and CO column 
over this temperature range may contribute to limiting 
the radial extent of the CO fundamental emission from 
the disk. 

The Spitzer HCN emission detected from T Tauri stars 
has been characterized as having a typical temperature 
of - 700K (Carr & Najita 2011; Salyk et al. 2011). If 
the emitting area is assumed to be the same as that of 
the water emission (i.e., ~ lAU), the line-of-sight col- 
umn density is typically ~ 10^^ cm~^ (Salyk et al. 2011). 



Fits made without that assumption yield smaller emit- 
ting areas (0.3 — 0.6 AU) and larger HCN columns (a few 
xlO^^ cm~^), although the lower columns of Salyk et al. 
are not ruled out (Carr & Najita 2011). In our models, 
HCN is abundant throughout the inner disk atmosphere, 
with a warm column of 10^'' — 10^^' cm^^ at 1 AU de- 
pending on the model parameters (Tables 3, 4, and 5), 
a range that encompasses the values inferred by Salyk 
et al. (2011). The larger line-of-sight columns preferred 
by Carr & Najita (2011) may not be too different from 
the warm (vertical) columns of our models if the HCN 
emission arises at disk radii closer to 0.5 AU, where the 
warm column is larger, and when disk inclination effects 
are taken into account. 

The similar temperatures found for the Spitzer HCN 
and water emission (Carr & Najita 2008, 2011; Salyk 
et al. 2008, 2011) suggest that the HCN emission arises 
from the same region of the disk atmosphere as the water 
emission. The smaller emitting area inferred for the HCN 
emission, compared to the water emission, in the best-fit 
models of Carr & Najita (2008, 2011) might seem odd 
in that context. The difference might indicate that the 
HCN emission is restricted to a smaller (possibly annu- 
lar) region than the water emission, perhaps because of 
the molecular abundance distribution or excitation con- 
siderations. In our models, the smooth variation of the 
warm HCN column with radius and the slow decline be- 
yond lAU (Figure 5) does not suggest an annular dis- 
tribution for the HCN. However, the radial range of the 
HCN emission may also be restricted by the high den- 
sities needed to excite the observed transitions and the 
decline in density with disk radius in the warm atmo- 
sphere (Figure 2). 

The C2H2 emission detected with Spitzer has been 
characterized as having a C2H2/HCN column density ra- 
tio of ~ 0.1 and a C2H2 column of lO^'"^ - lO^^cm-^ 
assuming that the C2H2 emission has the same tempera- 
ture and emitting area as the HCN emission (Carr & Na- 
jita 2008, 2011). Salyk et al. (2011) instead assumed that 
the C2H2 and HCN emission have the same emitting area 
as the water emission but allowed the C2H2 and HCN 
temperatures to differ; they inferred typical C2H2/HCN 
column ratios ~ 1 and columns of ^ 10^'' cm^^ for C2H2. 
In our reference model, the C2H2/HCN warm column ra- 
tio is ~ 0.01 and the C2H2 column density is ~ 10^^ cm~^ 
at 1 AU, both of which are smaller than indicated by 
the observations. However, for the case Ug ~ 7.07, the 
C2H2/HCN warm column ratio is ~ 0.1 and the C2H2 
column density is ^ lO^^cm"^ at lAU (Table 3), more 
similar to the reported values. 

The Spitzer CO2 emission has been fit with charac- 
teristic temperatures higher (Salyk et al. 2011) or lower 
than (Carr & Najita 2011) that of water. Carr & Na- 
jita (2011) find that the CO2 column density is not well 
constrained by the observations. Salyk et al. (2011) re- 
port CO2 column densities ^ 10^"'' — lO^'^cm^^, assum- 
ing the CO2 emission has the same emitting area as the 
H2O emission. In our models, the CO2 warm column is 
~ 10^5 - 10^*5 cm-2 at 1 AU (Tables 3, 4, and 5), similar 
to the reported values. 

As in GMN09, our models are unable to account for the 
OH emission from inner disks. The observed OH columns 
are 3 x lO^'^-Sx 10^^ cm'^ (Salyk et al. 2011; Carr & Na- 
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jita 2008), higher than in our models 10^''— 10^"'' cni^^; 
Tables 3, 4, and 5). The larger OH columns arc very 
likely the result of UV irradiation of the disk, which 
is not yet implemented here. UV irradiation can dis- 
sociate H2O to produce OH, potentially explaining the 
large OH columns that are observed (Bethell & Bergin 
2009; GMN09). Hot OH emission, which may resuh 
from the photodissociation of H2O (e.g., Najita et al. 
2010), has been reported from T Tauri disks (Carr & 
Najita 2011), consistent with this interpretation. Given 
the central role of the OH radical in the warm chem- 
istry of disk atmospheres (Figure 3), an enhanced OH 
abundance (e.g., from photodissociation) may lead to en- 
hanced abundances of species that are synthesized from 
OH. 

Our models address the question raised by Salyk et al. 
(2011) regarding the dominant carrier of nitrogen in the 
disk. Noting that the HCN abundance they infer does 
not account for all of the available nitrogen, Salyk et 
al. (2011) suggested NH3 or N2 as the possible carriers 
of the remaining nitrogen. They also pointed out how 
NH3 is unlikely to be the missing reservoir: given the 
upper limit on NH3 emission in the Spitzer spectra in 
the 10 lira region, they inferred that the column of warm 
NH3 is < lO^^cm"^, and N2 was suggested to be the 
dominant reservoir of nitrogen. Consistent with this in- 
terpretation, our model predicts warm NH3 columns of 
10^^ - 10^^ cm-2 at 1 AU (Tables 3, 4, and 5) and that 
most of the nitrogen is in N2 • 

4.1.2. Molecular Emission Trends 

Our models may also provide some context for the var- 
ious molecular emission trends that have been reported 
in the literature. In their study of a small sample of 
T Tauri stars, Carr & Najita (2011) found tentative evi- 
dence that the flux ratio of HCN/H2O emission increases 
with disk mass, increasing by a factor of ~ 4 over disk 
masses of 0.002 — 0.02 M©. They interpreted this result 
as a possible consequence of planetesimal or protoplanet 
formation. This interpretation arises from considering 
the transport of water in the solar nebula (Ciesla & Cuzzi 
2006). Because icy bodies ranging from micron, to me- 
ter, to kilometer and larger size migrate at different rates, 
they can transport water ice to the inner disk at greater 
or lesser rates compared to the gas, reducing or enhanc- 
ing the 0/C ratio of the inner disk and possibly affecting 
the chemistry of the inner disk. 

That is, when grains in the outer planet formation re- 
gion are small enough to couple well to the gas, water 
that is frozen on these grains accretes along with the gas 
into the inner planet formation region. The water is re- 
turned to the gas phase when the ices evaporate, and 
the 0/C ratio of the inner disk gas is representative of 
the disk as a whole. However, when the icy grains grow 
into meter-sized bodies, they migrate inward extremely 
rapidly relative to the gas. When they evaporate, they 
hydrate the inner disk, and in depositing water and oxy- 
gen potentially enhance the 0/C ratio there. In con- 
trast, if icy material in the outer planet formation region 
grows into planetesimal or larger bodies, such large bod- 
ies cease to migrate (and may go on to form planets). 
They thereby sequester water (and oxygen) in the outer 
planet formation region, leaving the inwardly migrat- 
ing material dehydrated and oxygen-poor (e.g., Ciesla 



& Cuzzi 2006). This situation might both reduce the 
water abundance in the inner disk and, by reducing the 
0/C ratio, favor the formation of hydrocarbons in the 
inner disk. 

We might therefore expect that early in the disk evolu- 
tion process (before significant grain growth has occured) 
the gaseous 0/C ratio of the inner planet formation re- 
gion will reflect the bulk 0/C ratio of the disk. At in- 
termediate times, the inner disk may have a higher 0/C 
than the disk as a whole. At late times, the inner disk 
may have a lower 0/C ratio than the disk as a whole. 
The timescale on which a disk transitions from "early" to 
"intermediate" to "late" depends on the disk mass: more 
massive disks would transition more rapidly through each 
phase. 

Thus, a coeval cluster of T Tauri stars with a range 
of disk masses may have inner disks with a range of gas- 
phase 0/C ratios, and we might expect the inner disks of 
T Tami stars to show a range of HCN/H2O flux ratios, as 
is observed (Carr & Najita 2011; Salyk et al. 2011). For 
T Tauri stars transitioning from the "intermediate" to 
"late" phase, HCN/H2O might increase with disk mass, 
possibly accounting for the tentative trend reported by 
Carr & Najita (2011). 

Our results support this scenario. As expected, a lower 
0/C ratio enhances the HCN abundance and decreases 
the H2O abundance in the inner disk. A very modest de- 
crease in the 0/C ratio, from 2.5 (our reference case) to 
~ 1.8 would raise the HCN/H2O ratio of warm columns 
by a factor of ~ 4. It would be interesting to measure 
other molecular ratios that might corroborate this pic- 
ture. We had earlier suggested the possibility of using 
C2H2/H2O as a diagnostic of the 0/C ratio (Section 3). 
Our results also suggest diagnostics of the "intermediate" 
phase of evolution. For the case of 0/C higher than the 
reference model (i.e., an 0/C ratio of 5), oxygen-bearing 
molecules such as H2O and O2 are enhanced and carbon- 
bearing molecules such as C2H2 and HCN are reduced in 
abundance relative to the reference case (Table 5). 

An alternative explanation for the tentative trend of 
HCN/H2O emission with disk mass reported by Carr & 
Najita (2011) may be an excitation effect, where HCN 
is better excited at the higher densities that would be 
present in a more massive disk. This scenario would 
require that the disk mass estimated from observations 
that probe large disk radii (> 20 AU from submillime- 
ter/millimeter observations) apply to the inner few AU 
of disks. These two possibilities (0/C ratio or excita- 
tion effect) could be tested by looking for trends with 
disk mass of emission from other molecules that are less 
sensitive to density. Detailed excitation modeling would 
also be useful in sorting out these possibilities. 

More perplexing is the report by Teske et al. (2010) of 
a possible trend of HCN flux with Lx, with HCN flux 
increasing on average by a factor of ^ 4 in the range 
Lx = 2 X 10^^ — 2 X lO'^^ergs"^. In our models, a higher 
Lx enhances the abundance of He"*" , which can break up 
CO and potentially enhance the HCN abundance. How- 
ever, our warm HCN columns vary negligibly with Lx in 
the T Tauri luminosity range (Lx = 10^^ — lO'"*^ ergs^^). 
It would be interesting to see if studies of larger samples 
of objects confirm a trend of HCN flux with Lx- 

There is a bigger difference between the warm HCN 
column for the reference case and for very low values of 
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Lx- For Lx = lO^^orgs"^, the warm HCN column is an 
order of magnitude larger than in the reference case (Ta- 
ble 4). As T Tauri stars are known to be strong emitters 
of X-rays at a much higher level than lO^^ergs"^, the 
ions produced by T Tauri X-ray emission appear to have 
the effect of limiting the HCN abundance and emission 
from disks. 

There is tentative evidence that HCN flux increases 
with stellar accretion rate, with HCN flux increasing 
from < 1 mJy-/zm to ~ 4mJy-/im over the range 10~^ — 
IQ-'^Mgyr-i (Teske et al. 2010). One possible explana- 
tion for this is that the increased HCN flux arises from 
increased heating. GMN09 had earlier discussed how en- 
hanced accretion-related mechanical heating in the disk 
atmosphere could enhance the water emission from the 
disk, primarily by increasing the column density of the 
warm atmosphere. We would expect a similar enhance- 
ment in HCN emission, since a warmer disk atmosphere 
would enhance the emission from other molecules as well 
as water, all other things being equal. 

Our results show that indeed increased accretion- 
related heating enhances the column of warm HCN in the 
disk atmosphere. The warm HCN column increases by a 
factor of - 3 (from 2.4 x IQi'^cm-^ to 7.7 x IQi'^cm-^) 
for values of = 0.01 to 1, with most of the increase 
in column occuring in the range a^, = 0.1 — 1. In the 
higher ah case, the warm HCN also occurs deeper in the 
atmosphere and at higher density where it may be better 
excited. If higher stellar accretion rates are accompanied 
by higher rates of mechanical heating in the disk atmo- 
sphere, both the increase in the warm column of HCN 
and the higher density region in which it is located may 
account for the increased HCN flux that is observed at 
higher stellar accretion rates. 

An alternative explanation is that the higher UV fluxes 
associated with accretion lead to greater grain photo- 
electric heating and enhanced warm molecular columns. 
Detailed modeling is needed to explore the impact of in- 
creased UV irradiation and whether the higher heating 
and photodissociation rates lead to higher or lower warm 
molecular columns. The result is likely to depend on the 
degree of grain settling in the atmosphere (e.g., Nomura 
et al. 2007). 

An alternative chemical explanation for the trend is 
that the higher UV fluxes dissociate N2 in favor of atomic 
N, which then enhances the formation and abundance 
of HCN. Such a scenario was discussed by Pascucci 
et al. (2009) who argued that the higher flux ratio of 
HCN/C2H2 in the spectra of T Tauri stars compared to 
those of brown dwarfs may result from the higher ac- 
cretion rates and consequently higher UV fluxes of T 
Tauri stars. In our models, a signiflcant fraction (ap- 
proximately half) of the nitrogen is already in atomic 
form in the warm atmosphere even in the absence of UV 
irradiation, with the remainder primarily in N2. Includ- 
ing the UV dissociation of N2 may then have a limited 
impact on the HCN abundance, although detailed calcu- 
lations are needed to be certain. 

4.1.3. Water Molecule Shielding and Heating 

The column density over which UV irradiation affects 
the disk atmosphere, either chemically (through pho- 
todissociation) or thermally (by depositing energy into 
the atmosphere), is governed by the abundance of UV- 



absorbing dust and molecules in the disk atmosphere. 
Molecules may provide the dominant source of UV opac- 
ity in T Tauri disk atmospheres given the high inferred 
rate of dust settling from disk atmospheres, as previously 
noted in the speciflc context of water by Bethell & Bergin 
(2009). 

Our results support this interpretation. Since water 
absorbs UV photons in the 1200-1800 A range with a 
cross-section of 10^^^ — lO"-"^^ cm^ (Parkinson & Yoshino 
2003), stellar UV photons encounter an optical depth 
of unity over a line-of-sight disk column density (to the 
star) of A^H20 ^ 10^^ — 10^*^ cm^^, or a vertical column 
density that is ~ 10 times smaller. The much larger total 
vertical water column densities that can be synthesized 
in our disk atmosphere models (» 10^* cm~^), the warm 
column component of which (~ 10^^ cm~^) is observed in 
emission from T Tauri disks (Carr k Najita 2008, 2011; 
Salyk et al. 2008, 2011), appear to be capable of shielding 
the lower disk atmosphere and the disk midplane from 
UV irradiation. 

As noted above, UV irradiation of the disk may also 
heat the disk surface through the absorption of UV pho- 
tons by water and other molecules. A simple estimate of 
the heating rate per unit volume due to the absorption 
of UV photons by water is 

Fuv = C*c- Euv XH2O n-H, 

where C* is the unshielded dissociation rate due to stel- 
lar UV, is the attenuation of the UV along a line 
of sight to a given position in the disk, and -Euv is the 
energy deposited as heat per dissociation. Laboratory 
experiments show that water is strongly dissociated in 
the wavelength range from 1200 I8OOA (e.g., Parkin- 
son & Yoshino 2003). Since the mean photon wave- 
length for dissociation is ^ 1500 A and the threshold is 
2340A, E'uv ^ 3eV. Bergin et al. (2003) showed that 
the stellar UV flux of the representative T Tauri star BP 
Tau at a distance of 100 AU is ~ 540 times the stan- 
dard Habing interstellar radiation field, or ^ 5.4 x 10^ 
larger at lAU. Multiplying the unattenuated interstel- 
lar photodissociation rate of 5.9 x 10~^°s~^ (Roberge 
et al. 1991) by 3 x 10^, after correcting for the use by 
Roberge et al. of the Draine interstellar radiation field, 
yields C* = 1-8 x lO'^s'^ 

Similarly, the X-ray heating rate is 

Fx = C,xe~^"" ExriH, 

where C,x is the unshielded X-ray ionization rate, Tx is 
the line of sight X-ray attenuation factor, and Ex is the 
heating per X-ray ionization. In our reference model, 
= 8 X 10~^ s~^ at the top of the atmosphere at 1 AU, 
and Ex 20 eV (GNI04). Substituting the given num- 
bers and assuming that .t(H20) = 10^*, the ratio of 
the UV water to X-ray heating rates is ^ 3e^'^"^/e^'^'=. 
This rough estimate confirms that UV water heating is 
potentially significant at the disk surface and that it is 
sensitive to the absorption of the UV by dust and by wa- 
ter itself. In the reference model of this paper at 1 AU, 
accretion heating in the neighborhood of the transition 
region dominates X-ray heating by a factor between 10 
and 100. A UV heating rate that is ^ 3 times the X-ray 
heating rate could make a quantitative difference, espe- 
cially where the transition occurs. 
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Thus water in the disk atmosphere may heat itself (cf. 
Bethell & Bergin 2009), a situation that would help to 
explain why the water emission observed with Spitzer 
has similar properties from source to source (Carr & Na- 
jita 2011). If water is a significant heat source at the 
disk surface, and if HCN is abundant at the disk surface 
(e.g., in the low Lx case or the case of ag = 7.07/im), 
UV heating by water may have a (thermal) role in ex- 
plaining a trcind bcitween HCN flux and stellar accretion 
rate. In addition to its potential role as a heat source, 
water is also known to be an important coolant based 
on its observed emission luminosity (Pontoppidan et al. 
2010). It would be important to explore the impact of 
water on both the heating and cooling of the disk atmo- 
sphere. Similar considerations apply to FUV absorption 
by H2, which is more abundant than H2O but absorbs 
over a more limited wavelength region. A detailed in- 
vestigation is needed to determine how and where these 
processes affect the disk properties. 

4.2. Outer Planet Formation Region 

The outer planet formation region, with its low dust 
temperature (e.g., < 150 K at 10 AU; Fig. 1), is more dif- 
ficult to treat with our current model becaiisc wc do not 
include freeze out and desorption processes. Moreover, 
the lower gas temperatures and densities in the outer 
planet formation region (Fig. 1) lead to slower chemical 
timescales. As a result, the steady state abundances that 
we calculate may be less valid, and transport processes 
(e.g., Willacy et al. 2006; Willacy & Woods 2009; No- 
mura et al. 2009; Heinzeller ct al. 2011), which arc not 
treated here, important. Nevertheless, some aspects of 
our results may bear on existing observations and future 
attempts to observe this region of the disk. 

Our simple models predict a dramatic decrease in the 
warm column of water beyond a few AU. At 10 AU and 
beyond, much of the water that is produced would freeze 
out at the low grain temperatures in the atmosphere. 
However, the model does predict a detectable water col- 
umn at smaller radii in the region of the disk atmosphere 
where the dust temperature is high enough that water 
would not freeze out. At 4 AU, the column of water above 
the point where the grain temperature reaches 150 K, 
125 K, and 100 K is 2 x 10^'^ cm.-'^, 2 x lO^^cm^^ and 
3 X 10^^ cm^^ respectively. Such a warm water column 
would likely produce detectable emission, in potential 
conflict with the Spitzer results, which find a small emit- 
ting area for the water emission of tt where r < 2 AU 
(Carr & Najita 2011; Salyk et al. 2011). This point was 
made by Meijerink ct al. (2009) based in part on the 
assumptions of the GMN09 model. 

Meijerink et al. suggested vertical transport as a way 
of resolving the potential conflict (see also Pontoppidan 
et al. 2010). In the scenario they proposed, water from 
the disk surface is transported, via vertical turbulent dif- 
fusion, further down in the disk to a region where the 
dust temperature is low enough that water can freeze 
out onto the grains. These icy grains settle to the disk 
midplane and are incorporated into large enough bodies 
that they are not rclofted to the disk surface. Meijerink 
et al. argue that this efl:ect (a "cold finger"), can account 
for the small emitting area inferred from the Spitzer wa- 
ter observations. They further suggested that the radial 
distribution of water emission at the disk surface might 



trace the location of the midplane snowline, which is of 
interest for planet formation. 

Since this picture, created to account for the proper- 
ties of the water emission observed with Spitzer, would 
also apply to other molecules with similar condensation 
temperatures, one way to test the picture is to see if such 
molecules have a similar radial distribution to that of wa- 
ter. One possible test molecule is HCN; its condensation 
temperature is similar to that of water (e.g., Walsh et al. 
2010) and our models predict it to remain abundant in 
the disk atmosphere out to 10 AU or more in the absence 
of freeze out. Thus far, the emitting areas inferred for 
HCN emission from the Spitzer data arc not sufficiently 
constrained to test this picture, as the HCN area can be 
smaller than, similar to, or greater than the water emit- 
ting area, given the uncertainties (Carr & Najita 2011). 
Searches for emission from cooler water and HCN, e.g., 
with the Herschel Space Telescope, may better address 
this issue. 

There are several considerations that would bear on 
whether our model results really motivate a "cold finger" 
for water or any other molecule, as differences between 
the model results and observations might be accounted 
for, to some extent, by some combination of the follow- 
ing considerations. Firstly, the disk atmosphere beyond 
a few AU may be cooler than we have inferred. If grains 
have settled less in the outer disk (e.g., than in the inner 
disk), we will have underestimated the gas-grain cooling 
there. Another concern is that we have neglected wa- 
ter cooling. Including such effects may reduce the warm 
water column and the molecular emission from the disk 
atmosphere without requiring vertical transport. Sec- 
ondly, the synthesis of molecules such as water may be 
less efficient than we have assumed. Some of the neu- 
tral reactions we have used have thermal barriers and 
their rate coefficients are very sensitive to temperature 
(e.g., eqs. 2 and 4) and may have been overestimated in 
the analyses of low temperature laboratory data. An- 
other important issue is the (uncertain) efficiency of H2 
formation on grains at high temperature, as discussed 
in Section 3. Thirdly, photodissociation may limit the 
amount of water that is present in the disk atmosphere. 
While processes such as photodissociation and a reduced 
efficiency of H2 formation on grains will affect both the 
inner and outer disks, they may have a greater relative 
impact on the outer disk (where the warm water col- 
umn is expected to be small even without these effects; 
Fig. 5) and may push the warm water column below an 
observable level in the outer disk. 

We can test the first point above by looking for molec- 
ular emission from molecules other than water. Our 
model predicts a thermal-chemical transition out to be- 
yond 40 AU and that CO is abundant in the disk at- 
mosphere over that range of radii. The presence or ab- 
sence of any resulting warm CO emission (e.g., with the 
Herschel Space Telescope) would probe whether or not a 
warm disk atmosphere is present at radii beyond a few 
AU. 

How well can the models presented here really identify 
chemical signatures of planetcsimal formation, one of the 
longer-term motivations for our study? Our variation 
of the 0/C ratio in the inner planet formation region 
of the disk may capture some of the relevant behavior 
of the scenario discussed in section 4.1.2. However, a 
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detailed freeze out and desorption model and the inclu- 
sion of transport processes (e.g., as explored by other au- 
thors) is likely needed to predict which species remain in 
the disk atmosphere and at what abundance. There are 
a number of uncertainties in constructing such a model 
(see e.g., Walsh et al. 2010 for a discussion of desorption 
processes). This is an important area for future work. 

One tracer of the disk atmosphere that is independent 
of many of the above concerns is HI 21 cm radiation, 
which is of interest because one feature of X-ray irradi- 
ated disk atmospheres is a significant HI component. Our 
models indicate an HI column of 2 x 10^^ — 7 x 10^^ cm~^ 
in the 1-20 AU region from each face of the disk. A simi- 
lar column extending to 100 AU would imply an HI mass 
of ~ 10^* gm, a value that is consistent with current HI 
upper limits but may be detectable with future observa- 
tional facilities (Kamp ct al. 2007). 

5. SUMMARY AND FUTURE DIRECTIONS 

As discussed in the introduction, one goal of our mod- 
eling effort was to explore the potential role of pro- 
cesses such as grain settling, X-ray irradiation, accretion- 
related mechanical heating, and the 0/C ratio in de- 
termining the thermal-chemical properties of disk atmo- 
spheres. Another was to explore how well a model that 
excludes UV irradiation can do in accounting for the ex- 
isting observations of inner disk atmospheres. We also 
aimed to make detailed comparisons with the observa- 
tions in order to determine the changes needed to im- 
prove our model. 

To address these goals, we used an expanded version of 
our thermal-chemical model of disk atmospheres to ex- 
plore the sensitivity of the properties of the atmosphere 
to model parameters. We find that the warm columns 
(200-2000 K) of many molecular species are sensitive to 
grain settling and the efficiency of accretion heating, pri- 
marily because these parameters affect the depth of the 
warm molecular layer at the disk surface. We also showed 
how certain model parameters (grain growth, stellar X- 
ray luminosity, mechanical heating, 0/C ratio) affect the 
molecular abundances of disk atmospheres, and which 
species are sensitive to which parameters. 

We find many areas of agreement with the observa- 
tions. The model parameter variations we considered can 
account for many of the properties of molecular species 
that have been detected with Spitzer and ground-based 
near-infrared spectroscopy (i.e., their typical tempera- 
tures, warm columns, and emitting areas). The sensitiv- 
ity of warm columns to parameters such as grain settling 
and accretion heating may account, at least in part, for 
the large range in the molecular emission fiuxes that have 
been observed from T Tauri disks. 

The dependence of the warm columns on model param- 
eters such as grain growth, accretion heating, and 0/C 
ratio, may help to explain the trends reported in the lit- 
erature between the detection rate of molecular emission 
and mid-infrared color (Salyk et al. 2011) and molecu- 
lar emission fiux with stellar accretion rate (Teske et al. 
2011) and disk mass. Regarding the last, we find that 
the ratio of HCN/H2O increases with a decreasing 0/C 
ratio; a decreasing 0/C ratio in the inner disk is a po- 
tential effect of planetesimal and protoplanet formation 
in the outer disk, which likely occurs at an accelerated 
rate in higher mass disks (Carr & Najita 2011). 



While the parameter variations we explored can ac- 
count for the average values and ranges of the observed 
molecular emission properties of disks, a future challenge 
will be to explain the properties of individual sources 
with a single set of model parameters. Our model suffers 
from some obvious limitations; e.g., our low OH abun- 
dance points out the need for including UV irradiation 
and photodissociation in our model. The results of our 
steady-state chemical model are also likely less valid in 
the outer planet formation region of the disk because 
of the longer chemical timescales there and because we 
have not included freezeout, desorption and transport 
processes. We aim to address these challenges in future 
work. 
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Table 1. Elemental Abundances 



Element 


Abundance 


Depletion 


H 


1.00 


1.0 


He 


0.10 


1.0 


C 


1.40 X 10-"* 


2.0 


N 


6.00 X 10-5 


1.1 





3.50 X 10-^ 


1.4 


Ne 


6.90 X 10"'^ 


1.2 


Na 


2.31 X 10"^ 


7.5 


Mg 


1.00 X 10"^ 


40 


K 


8.57 X 10-9 


12 


Si 


1.68 X 10"*^ 


19 


S 


1.40 X 10-5 


1.0 


Ar 


1.51 X 10-6 


1.7 


Fe 


1.75 X 10-"^ 


180 



Table 3. Warm Columns (200-2000 K)t 



Oih 


1.0* 


1.0 


0.01 


Ug 


0.707* 


7.07 


0.707 


H2 


3.91E+21 


2.40E+22 


1.33E+21 


H+ 


3.20E+12 


3.G0E+10 


9.56E+15 


H+ 




1.92E+13 


1.45E+12 


2.78E+13 


He+ 


7.75E+12 


5.00E+11 


2.28E+14 





1.70E+17 


3.70E+16 


1.76E+17 


OH 


4.57E+14 


6.83E+13 


1.31E+15 


H2O 


5.40E+17 


5.29E+18 


1.96E+17 


H3O+ 


4.22E+13 


1.36E+13 


7.50E+13 


O2 


3.72E+17 


1.73E+18 


1.05E+17 


CO 


1.09E+18 


6.55E+18 


4.30E+17 


C02 


1.26E+15 


l.lOE+16 


2.25E+14 


HCO 


1.33E+16 


8.53E+16 


3.54E+14 


HCO+ 


2.28E+13 


1.49E+12 


4.60E+13 


C 


2.68E+14 


3.33E+13 


9.34E+15 


CH3 


3.78E+09 


2.11E+10 


1.41E+09 


CH4 


9.59E+11 


1.28E+14 


7.35E+10 


C2H2 


1.12E+13 


5.53E+14 


1.36E+12 


c+ 


4.11E+12 


1.33E+11 


1.22E+15 


CH^ 


5.48E+11 


1.69E+11 


1.41E+11 


C2Ht 


3.79E+08 


3.08E+08 


3.73E+07 


N 


9.01E+16 


1.24E+17 


4.59E+16 


N2 


1.92E+17 


1.37E+18 


7.10E+16 


NH3 


3.42E+15 


1.14E+16 


8.40E+14 


HON 


8.05E+14 


7.29E+15 


2.50E+14 


NH;j" 


5.30E+11 


3.50E+10 


7.44E+11 


N2H+ 


1.53E+12 


1.42E+11 


2.73E+12 


HCNH+ 


1.96E+11 


3.67E+11 


3.20E+11 


NO 


7.29E+14 


1.23E+14 


1.55E+15 


S 


1.28E+16 


8.65E+16 


1.42E+15 


s+ 


2.42E+13 


9.25E+11 


6.04E+14 


so 


7.23E+14 


3.78E+15 


3.75E+14 


S02 


9.79E+16 


5.70E+17 


4.21E+16 


cs 


l.OlE+14 


3.87E+14 


1.04E+13 


^ Gg in fim and warm 


columns in 


cm ^. 



* Reference case. 



Table 2. Model Parameters 



Parameter 


Symbol 


Ref Value 


Variation Range 


stellar Mass 


M* 


0.5Mo 




Stellar Radius 


R* 


2Rq 




Stellar Temperature 




4000 K 




Disk Mass 


Md 


0.005Mo 




Disk Accretion Rate 


M 


lO-*M0yr-i 




Grain Size 


ag 


0.707/im 


0.707 - im^im 


X-ray Luminosity 


Lx 


2 X 103"ergs"i 


1026 _ io3iergs"^ 


X-ray Temperature 


Tx 


IkeV 


1-5 keV 


Mechanical Heating 


ah 


1.0 


0.01 1.0 


0/C Abundance Ratio 


0/C 


2.5 


0.25 - 5 


Radial Distance 


r 


lAU 


0.25-40 AU 



16 



Table 4. Dependence of Warm Columns (200-2000 K) on 





1.0c26 


1.0e27 


1.0c28 


1.0e29 


1.0e30 


2.0e30* 


1.0e31 


Ho 




Szip,-I_91 





«zlF-l-91 


Q 



SzlF-l-91 


Q 


7QF-I-91 


"3 



R7F-I-91 


Q 



U -LLJ-/T^Z ± 


Q 



q«F-l-91 
y OL^n^z ± 


H+ 


9 




Q 


1 7F,-l-nS 


■3 



7nF4-nQ 

1 wj_;^wy 


A 


'17F-4-1 n 


a 




49F-I-1 1 


9 
z 


■^DF-I-I 9 

OULl/T^ -L Z 


9 
z 


R«F-I-1 

UOL-^n^ ±0 


H+ 

3 


1 

J_ 


9fiP-l-flQ 


1 
1 


zL'^F-l-1 n 


1 

X 




1 


■^OF^I 9 


y 


ZL9F-I-1 9 

^Zllin^ -LZ 


1 

-L 


fi9F^1 

UZilL/^ 1 


zl 

4: 




He+ 


A 




A 

4: 


S'^F-l-OQ 
oorijn^uy 


c; 


1 SF^I 


zl 


^1 F ^ 1 1 





•^flF-H 9 

t-)U-LJn^ -LZ 


f! 
u 


ynp 1 1 9 
oy-L-/^ 1 z 


9 
z 


U4;-LJn^ -LO 





u 




u 




f! 
u 


"^1 F+1 ^ 


Zl 

4: 




1 
± 


1 7F-I-1 7 


1 

± 


4nF+1 7 


9 
z 


^7F-I-1 7 


OH 






9 


nqF-i-1 9 


1 

X 




R 
U 




9 
z 




Q 


yU-EJT^ -L^ 


1 

-L 




H2O 


q 


zIc:F,-I-1 7 


a 




i^nF-H 7 


7 


fi7F,-l-1 7 


R 
u 


4'?F,-I-1 7 


K 



49F,-I-1 7 


c: 
(J 




R 
U 


9c;f,-I-1 7 

Z(J -LJn^ ± / 


H3O+ 







Q 





1 
± 




u 


1 ^F+1 9 


9 
z 




Q 





1 
± 


99F-I-M 


O2 


9 


t^QF-Ul 7 


9 


snF4-i 7 


Q 


1 7F-I-1 7 


Q 



'i9F4-1 7 


Q 


4SF-I-1 7 


Q 






ztztF-l-1 7 


CO 


1 

J. 


nnF-n 8 

UWJ_;T^ ±0 


1 

± 




1 


n^iF-i-i s 


1 

± 


U'il-JT^ -LO 


1 

± 


UZ-LL/^ -LO 


1 

-L 


U-LL-J^ ±0 


1 

± 


±U-i--/n^ -LO 


C02 


I 


77F,+1 ^1 

1 J -Li 1 -1- ' J 


\ 


4QF+1 




'?'?F-)-1 ^ 

O0-Li~ J- 




J- rJ-LJ 1 -L U 






I 


99F+1 ^ 

Z Z J_J 1^ _L 


X 


fi3F+1 ^ 

UOJ-J 1^ XfJ 


HCO 







Q 
y 




q 
y 




q 
y 




1 
± 


^UXJn^ -LU 


1 
± 


'^OF-i-l f\ 

OUJ-iJT^lU 


1 

X 




HCO+ 


Q 


OUJ— ;T^UO 


1 

± 


4ftF,_Ll 


9 




1 
± 


7QF-I-1 9 


1 

± 




1 
± 








c 


Q 




A 

4: 


±(jj_jn^ -1- -1- 


9 




1 

± 




1 
± 


':i4F-kl4 


9 

z 


QDF^l 4 

OU-LJ 1 -L4: 


7 
1 


(^4F-I-14 


CH3 


1 

J. 




7 




1 

± 




A 









Q 
u 




c; 




CH4 


u 




q 

y 




9 


1 9F,-I-1 


(J 


9c;F-I-1 9 


1 

± 




1 

± 


nOF-l-l 9 


c; 


4SF-H 1 

tOUJn^ -L -L 


C2H2 




7QF+1 fi 


4 


O'^F+14 

VJ J-J -1- t: 


a 


fifiF+1 3 


2 


1 flF+l 3 


Q 


70F+1 9 


I 


nfiF+1 3 


X 


77F-I-1 

1 1 J_jn^ -LO 


c+ 







1 
± 




1 


97F4-1 


1 
± 




1 
1 




Q 


•^1 F-Ul 9 
01 J_/n^ 1 z 


1 
± 


fi9F-Ul 

U^XJn^ -LO 


CH^ 


8 






R 
u 




9 




q 

y 




Q 





4 


7fiF-l-1 1 


1 

± 


OAF-H 9 

utn/n^ -LZi 


C2H+ 


1 
± 




9 









fi 
u 


n7F^n7 


1 

± 




Q 



Q OF -I- OR 


u 


7QF-kns 

( O-Ljn^wo 


N 


9 








1 




Q 





7 
1 


91 F-k1 B 

Zl J. J_J -L U 


a 



91 F-l-1 


1 

± 


14F_k1 7 

-L '4t.l_J -L 1 


No 


1 

± 


Q7P_I_1 7 


9 


91 F-Ul 7 


9 


99F-I-1 7 


9 


noF-i-i 7 


1 

L 


R4F-I-1 7 


1 
± 


77F-I-1 7 


1 

X 


£70Ji/n^i / 







Q7P_Ll R 

( rj-j- 10 


7 




1 

1 




9 
z 


^7F^1 ^ 
i Hz-f- 10 






Q 



UOrj-r 10 


Q 


1 9F-U1 ^ 


HCN 


fl 











1 

J. 


97F-I-1 ^ 


q 
y 


1 1 F-l- 1 4 


7 




7 


(^QF-l-l 4 


a 



49F-k14 


NH^ 


1 


08E+08 


7 


65E+08 


3 


77E+09 


3 


14E+10 


2 


63E+11 


4 


50E+11 


1 


34E+12 


N2H+ 


7 


34E+07 


1 


34E+09 


1 


65E+10 


1 


36E+11 


7 


97E+11 


1 


30E+12 


3 


72E+12 


HCNH+ 


4 


55E+10 


1 


04E+11 


6 


30E+10 


6 


68E+10 


1 


30E+11 


1 


74E+11 


5 


51E+11 


NO 


4 


99E+11 


3 


56E+12 


2 


72E+13 


1 


07E+14 


4 


13E+14 


6 


32E-I-14 


1 


78E+15 


S 


1 


06E+15 


3 


48E+15 


8 


OlE+15 


1 


30E+16 


1 


25E+16 


1 


25E+16 


1 


50E+16 


s+ 


1 


44E+08 


1 


39E+09 


4 


21E+10 


9 


53E+11 


1 


06E+13 


2 


OlE+13 


9 


82E+13 


so 


4 


79E+13 


2 


18E+14 


3 


76E+14 


4 


53E+14 


5 


80E+14 


6 


79E+14 


1 


33E+15 


S02 


8 


99E+16 


9 


15E+16 


9 


18E-I-16 


9 


17E-I-16 


9 


06E+16 


8 


96E-I-16 


1 


02E+17 


cs 


6 


98E+14 


6 


98E+13 


4 


66E+13 


5 


60E+13 


8 


39E+13 


9 


64E+13 


1 


51E+14 



^ Lx in ergs ^ and warm columns in cm 
* Reference case. 
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Table 5. Dependence of Warm Columns (200-2000K) on 0/C Ratipt 



0/C 


0.25 


0.50 


1.00 


1.50 


2.00 


2.5* 


5.00 


H2 


3 


85E+21 


3.88E+21 


3 


88E+21 


3 


91E+21 


3 


91E+21 


3 


91E+21 


3 


93E+21 


H+ 


9 


86E+13 


4.40E+14 


9 


21E+14 


1 


09E+13 


4 


86E+12 


3 


20E+12 


1 


44E+12 




3 


40E+13 


4.06E+13 


5 


08E+13 


3 


41E+13 


2 


41E+13 


1 


92E+13 


1 


llE+13 


He+ 


2 


18E+13 


2.95E+13 


1 


42E+13 


1 


n5E+13 


8 


81E+12 


7 


75E+12 


5 


64E+12 





1 


61E+12 


4.85E+12 


1 


17E+17 


1 


66E+17 


1 


76E+17 


1 


70E+17 


1 


45E+17 


OH 


2 


78E+11 


8.09E+11 


4 


31E+13 


2 


67E+14 


3 


93E+14 


4 


57E+14 


5 


82E+14 


H2O 


3 


68E+10 


1.64E+11 


9 


15E+14 


1 


48E+17 


3 


45E+17 


5 


40E+17 


1 


62E+18 


H3O+ 


1 


26E+08 


6.45E+08 


2 


28E+12 


2 


78E+13 


3 


68E+13 


4 


22E+13 


6 


24E+13 


O2 


2 


2lE+n8 


1.50E+09 


5 


74E+13 


2 


42E+16 


1 


95E+17 


3 


72E+17 


1 


25E+18 


CO 


1 


79E+16 


3.69E+16 


8 


78E+17 


1 


07E+18 


1 


n8E+18 


1 


nqE+18 


1 


llE+18 


C02 


3 


65E+12 


1.60E+13 


5 


85E+14 


1 


27E+15 


1 


27E+15 


1 


26E+15 


1 


20E+15 


HCO 


2 


54E+17 


5.12E+17 


7 


91E+16 


3 


17E+16 


1 


94E+16 


1 


33E+16 


3 


70E+15 


HCO+ 


5 


40E+12 


1.39E+13 


7 


16E+13 


5 


OOE+13 


2 


97E+13 


2 


28E+13 


1 


47E+13 


c 


5 


45E+17 


4.25E+17 


1 


12E+17 


7 


23E+14 


3 


83E+14 


2 


68E+14 


9 


OlE+13 


CH3 


2 


36E+11 


9.99E+10 


2 


37E+10 


3 


«qE+09 


3 


SOE-i-nq 


3 


78E+09 


6 


5qE+09 


CH4 


1 


88E+11 


1.39E+11 


3 


03E+11 


1 


2'^E-l-l 2 


1 


n3E-|-12 


9 


59E+11 


7 


56E+11 


C2H2 


1 


24E+17 


5.64E+16 


9 


34E+15 


2 


q2E+1 3 


1 


53E+1 3 




12E+13 


3 


t If 1 1 J 1 _L 


c+ 


1 


80E+13 


1.53E+13 


8 


14E+14 


3 


77E+13 


7 


iiE+12 


4 


llE+12 


1 


16E+12 


CH+ 


2 


53EJ.12 


2.33E+12 


1 


nnE+12 

\J\J l—J \ _L ^ 


9 


03E+11 


6 


84E+11 


5 


48E+11 


2 


91E+11 


C2H+ 


1 


51E+12 


8.55E+11 


3 


77E+10 


1 


53E4-09 


6 


8nE+08 


3 


7qp 1 08 


5 


23E+07 


N 


3 


37E+17 


3.76E+17 


3 


19E+17 


1 


57E+17 


1 


lOE+17 


9 


OlE+16 


5 


34E+16 


N2 


4 


90E+16 


3.48E+16 


7 


32E+16 


1 


59E+17 


1 


82E+17 


1 


92E+17 


2 


12E+17 


NH3 


4 


76E+15 


5.67E+15 


1 


05E+14 


1 


20E+15 


2 


87E+15 


3 


42E+15 


3 


fiqE-i-15 


HCN 


1 


93E+16 


9.55E+15 


8 


13E+15 


2 


lOE+15 


1 


17E+15 


8 


05E+14 


2 


90E+14 


NH+ 


1 


17E+12 


1.45E+12 


1 


96E+12 


1 


28E+12 


7 


56E+11 


5 


30E+11 


2 


35E+11 


N2H+ 


4 


65E+11 


3.87E+11 


4 


65E+11 


2 


05E+12 


1 


81E+12 


1 


53E+12 


9 


69E+11 


NH+ 


1 


17E+12 


1.45E+12 


1 


96E+12 


1 


28E+12 


7 


56E+11 


5 


30E+11 


2 


35E+11 


HCNH+ 


1 


08E+13 


9.80E+12 


2 


96E+12 


4 


42E+11 


2 


60E+11 


1 


96E+11 


1 


23E+11 


NO 


8 


48E+10 


3.24E+11 


3 


89E+13 


4 


25E+14 


6 


28E+14 


7 


29E+14 


1 


05E+15 


S 


1 


OlE+17 


1.03E+17 


9 


06E+16 


1 


51E+16 


1 


36E+16 


1 


28E+16 


1 


16E+16 


s+ 


5 


90E+15 


6.06E+15 


9 


44E+15 


6 


26E+13 


3 


04E+13 


2 


42E+13 


1 


83E+13 


so 


5 


94E+10 


2.35E+11 


9 


94E+13 


4 


32E+14 


5 


78E+14 


7 


23E+14 


1 


31E+15 


S02 


3 


54E+07 


3.96E+08 


1 


04E+16 


9 


59E+16 


9 


72E+16 


9 


79E+16 


9 


94E+16 


cs 


2 


55E+15 


1.44E+15 


1 


55E+14 


9 


78E+13 


1 


07E+14 


1 


OlE+14 


5 


20E+13 



^ Warm columns in cm 



* Reference case. 
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Table 6. Dependence of Warm Columns (200-2000K) on Disk Radius^ 



R (AU) 


0.25 


0.5 


1 


* 


2 




4 




10 


H2 


4 


85E+23 


1 


54E+22 


3 


91E+21 


2 


35E+21 


1 


68E+21 


7 


76E+20 


H+ 


9 


89E+12 


4 


81E+12 


3 


20E+12 


2 


29E+12 


1 


74E+12 


2 


40E+12 





3 


\_/tJJ — ' 1 J-fJ 


2 


R3E+13 


1 


92E+13 


1 


49E+13 


1 


16E+13 


9 


47E+12 





1 


59E+16 


4 


45E+17 


1 


70E+17 


9 


61E+16 


7 


29E+16 


7 


95E+16 


OH 


8 


OlE+14 


8 


lOE+14 


4 


57E+14 


3 


34E+14 


2 


94E+14 


4 


20E+14 


H2O 


1 


66E+20 


9 


75E+17 


5 


40E+17 


4 


17E+17 


3 


OOE+17 


5 


40E+16 


H3O+ 


1 


30E+14 


6 


21E+13 


4 


22E+13 


3 


20E+13 


2 


26E+13 


6 


92E+12 


O2 


6 


54E+18 


2 


08E+18 


3 


72E+17 


1 


85E+17 


1 


28E+17 


7 


81E+16 


CO 


1 


3nE+20 


4 


28E+18 


1 


09E+18 


6 


56E+17 


4 


74E+17 


2 


25E+17 


C02 


1 


3nE+16 

fJKJJ — 1 1 J. \J 


1 


77E+15 


1 


26E+15 


1 


09E+15 


7 


73E+14 


2 


22E+13 


HCO 


3 


21E+16 


1 


72E+16 


1 


33E+16 


1 


llE+16 


5 


62E+15 


9 


52E+11 


HCO+ 


1 


«qE+13 


3 


42E+13 


2 


28E+13 


1 


67E+13 


1 


41E+13 


1 


82E+13 


C 


6 


60E+14 


3 


61E+14 
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Fig. 1. — X-ray ionization rate (top panel), gas and dust temper- 
ature (middle panel), and abundances of atomic hydrogen, molecu- 
lar hydrogen, and electrons (bottom panel) as a function of vertical 
column density and radius in the reference model. 
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Fig. 2. — (a) Abundances of various species in the reference model 
at 0.25 AU (left), lAU (middle), and 2AU (right). 
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Fig. 2. — (b) Abundances of various species in the reference model 
at 4AU (left), 10 AU (middle), and 20 AU (right). 




Fig. 3. - Schematic chemical network showing how fast neutral 
reactions of OH with C, O, and N atoms generate oxides that 
provide pathways to CO2, O2, and N2. 
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Fig. 4. — Gas (solid) and dust (dashed) temperatures (top panel), 
electron (dashed) and H2 (solid) abundances (middle panel), and 
CO (dashed) and water (solid) abundances (bottom panel) in the 
reference case (black), and for = 1, Qg = 7.07 (red) and = 
0.01, ag = 0.707 (blue). 
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Fig. 5. — Warm columns of various species as a function of disk 
radius in the reference case. 



